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ABSTRACT 

The  dynamics  of  0(  P)  +  CO  collisions  at  a  hyperthermal  collision  energy  of  83  kcal  mol  have 
been  studied  with  a  crossed  molecular  beams  experiment  and  with  quasiclassical  trajectory 
calculations  on  computed  potential  energy  surfaces.  In  the  experiment,  a  rotatable  mass 
spectrometer  detector  was  used  to  monitor  inelastically  and  reactively  scattered  products  as  a 
function  of  velocity  and  scattering  angle.  From  these  data,  center-of-mass  (c.m.)  translational 
energy  and  angular  distributions  were  derived  for  the  inelastic  and  reactive  channels. 
Isotopically  labeled  C180  was  used  to  distinguish  the  reactive  channel  (160  +  C180  — >  16OC  + 
l80)  from  the  inelastic  channel  (160  +  C180  — >  160  +  C180).  The  reactive  16OC  molecules 
scattered  predominantly  in  the  forward  direction — i.e.,  in  the  same  direction  as  the  velocity 
vector  of  the  reagent  O  atoms  in  the  c.m.  frame.  The  c.m.  translational  energy  distribution  of  the 
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reactively  scattered  160C  and  180  was  very  broad,  indicating  that  16OC  is  fonned  with  a  wide 
range  of  internal  energies,  with  an  average  internal  excitation  of  60  percent  of  the  available 
energy.  The  c.rn.  translational  energy  distribution  of  the  inelastically  scattered  C  O  and  O 
products  indicated  that  an  average  of  15  percent  of  the  collision  energy  went  into  internal 
excitation  of  C180,  although  a  small  fraction  of  the  collisions  transferred  nearly  all  the  collision 

i  o 

energy  into  internal  excitation  of  C  O.  The  theoretical  calculations,  which  extend  previously 
published  results  on  this  system,  predict  c.m.  translational  and  angular  distributions  that  are  in 
near  quantitative  agreement  with  the  experimentally  derived  distributions.  The  theoretical 
calculations,  thus  validated  by  the  experimental  results,  have  been  used  to  derive  internal  state 
distributions  of  scattered  CO  products  and  to  probe  in  detail  the  interactions  that  lead  to  the 
observed  dynamical  behavior. 


I.  INTRODUCTION 

Hyperthennal  oxygen-atom  collisions  are  abundant  on  and  around  space  vehicles  that 
travel  at  high  velocities  through  the  rarefied  atmosphere  of  the  Earth  at  altitudes  of  150-700  km. 
In  this  region  of  the  upper  atmosphere,  atomic  oxygen  (and,  with  a  lower  flux,  molecular 
nitrogen)  can  collide  with  space  vehicles  and  their  exhaust  streams  at  relative  velocities  in  the 
vicinity  of  8  km  s-1.1"4  These  high  relative  velocities  may  lead  to  gas-phase  collisions  with  more 
than  100  kcal  mol  1  of  collision  energy  in  the  center-of-mass  (c.m.)  reference  frame.4  Such 
collisions  may  result  in  large  energy  transfers  into  product  internal  degrees  of  freedom  and  the 
opening  of  reaction  channels  with  high  barriers. 
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Hyperthermal  collisions  between  O(  'P)  and  the  common  combustion  product,  CO(  'l  ), 
have  received  a  great  deal  of  attention.  When  the  reactants  and  products  are  in  their  electronic 
ground  states,  the  possible  outcomes  of  such  collisions  are: 

O  +  CO(v,j)  — >  O  +  CO(v',j')  inelastic  AH°o=0  (1) 

O  +  CO(v,j)  — »■ OC(v',j')  +  O  reactive  AH°o=0  (2) 

O  +  CO(v,j)  — ►  C  +  02(v',j')  reactive  AZ7°0  =  256  kcal  mol  1  (3) 

At  a  typical  relative  velocity  of  8  km  s-1,  the  c.m.  collision  energy  of  O  +  CO  is  77.9  kcal  mol-1. 
If  exhaust  gases  were  directed  at  3  km  s  1  into  the  ram  direction  of  a  spacecraft,  the  relative 
velocity  of  O  and  CO  could  be  ~1 1  km  s_1,  in  which  case  the  c.m.  collision  energy  would  be  147 
kcal  mol  1 .  Achieving  the  collision  energy  needed  for  the  formation  of  molecular  oxygen  (3) 
would  require  relative  velocities  that  are  unattainable  under  any  reasonable  operating  conditions. 
Therefore,  the  only  channels  that  are  accessible  at  the  collision  energies  relevant  to  exhaust  gases 
from  space  vehicles  in  the  upper  atmosphere  are  (1)  and  (2).  Upschulte  and  Caledonia5 
measured  infrared  emission  of  vibrationally  excited  CO  following  collisions  between  atomic 
oxygen  and  CO  at  relative  velocities  of  approximately  8  km  s_1,  and  they  reported  a  cross  section 
for  vibrational  excitation,  ^/cr(v  =  /) ,  for  CO(v)  of  7.3  x  1CT17  cm2.  From  their  measured 

i= 1 

spectra,  they  concluded  that  the  CO  molecules  were  highly  vibrationally  excited,  up  to 
CO(v=ll).  Their  result  presumably  contained  contributions  from  both  inelastic  and  reactive 
scattering.  A  space-based  study  by  Green  and  coworkers,6  where  optical  emission  spectra  were 
measured  from  a  remote  spacecraft  after  CO  gas  was  released  from  a  canister  in  low  Earth  orbit, 
determined  that  bright  infrared  emissions  in  both  fundamental  and  overtone  vibrational  bands 
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resulted  from  the  highly  energetic  interaction  of  the  released  CO  with  ambient  O('P)  at  relative 
velocities  of  ~8  km  s_1.  Their  measured  spectra  showed  substantial  rotational  and  vibrational 
excitation  in  the  CO  product,  closely  matching  the  results  of  Upschulte  and  Caledonia.  Green  et 
al  6  estimated  the  cross  section  to  create  vibrationally  excited  CO(v)  to  be  on  the  order  of  1 0  1 6 
cm2,  in  reasonable  agreement  with  the  results  of  Upschulte  and  Caledonia. 

Braunstein  and  Duff7  computed  the  potential  energy  surfaces  of  the  three  lowest 
electronic  triplet  states  of  CO?  (one  A’  state  and  two  A"  states),  which  correlate  with  Of  P)  + 
CO('X+).  In  addition,  they  ran  classical  trajectories  on  these  surfaces.  Fig.  1  shows  the 
calculated  potential  energy  curves  of  the  ground  singlet  state  of  CO?^')  and  the  three  lowest 
triplet  states,  with  the  O-C-O  angle  equal  to  120°.  The  energy  shown  here  is  a  function  of  one 
of  the  CO  bond  lengths,  while  the  length  of  the  other  CO  bond  is  fixed  at  1 .2  A.  For  the  lowest 
triplet  state  at  its  transition  state  geometry,  which  is  slightly  different  than  shown  in  the  figure, 
the  barrier  to  O-atom  exchange,  Eq.  (2),  was  calculated  to  be  4.6  kcal  mol  '.  Although  coupling 
to  the  singlet  ground  state  is  spin-forbidden  from  the  three  lowest-lying  triplet  states,  the  singlet 
potential  energy  surface  (the  ground  state  of  C02,  X  *Eg+)  crosses  the  triplet  surface  and  leads  to 
a  higher  dissociation  limit,  with  products,  O('D)  +  CO(  'l+),  about  46  kcal  mol  1  above  the 
dissociation  limit  that  leads  to  0(  P)  +  CO(  E  ).  Intersystem  crossing  (ISC)  is  thus  possible, 
although  spin  conservation  is  generally  favored  for  systems  containing  first-row  atoms,  and  other 

O  Q 

studies  involving  hyperthermal  O-atom  reactions  showed  no  evidence  for  ISC.  ’  Quasiclassical 
trajectory  calculations  on  the  computed  triplet  surfaces7  agreed  well  with  available  kinetic  data 
on  vibrational  relaxation  of  CO(v=l  — >  v=0)  by  0(  P),  except  at  low  temperature.  The  shape  of 
the  emission  spectra  of  CO  based  on  the  calculations  for  collisions  at  a  relative  velocity  of  8  km 
s“‘  were  in  fair  agreement  with  the  experimental  observations  of  Upschulte  and  Caledonia  and 
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with  the  spaced-based  observations  of  Green  et  al.  However,  the  magnitude  of  the  calculated 
cross  section  for  the  formation  of  vibrationally  excited  CO  was  about  an  order  of  magnitude 
above  the  two  experimental  values.  A  subsequent  modeling  study10  suggested  that  the  difference 
in  cross  sections  between  the  measured  Upschulte  and  Caledonia  result  and  the  calculations  of 
Braunstein  and  Duff  was  the  result  of  multiple-collision  effects  not  fully  accounted  for  in  the 
measurements,  although  this  suggestion  has  not  been  experimentally  confirmed.  The  calculated 
trajectories  were  studied  to  leam  mechanistic  details  of  the  collisions  that  lead  to  vibrational 
excitation  through  inelastic  and  reactive  collisions.  It  was  found  that  these  processes  occur  when 
the  reagent  O  atom  approaches  the  carbon  end  of  CO  on  the  two  low-lying  triplet  surfaces,  A' 

-5 

and  1  A",  which  have  small  wells,  at  O-C-O  angles  between  80°  and  140°. 

We  describe  here  a  combined  theoretical  and  crossed-molecular-beams  study  of  0(  P)  + 
CO  collisions  at  a  hyperthennal  collision  energy  near  80  kcal  mol-1  (Eco\\  =  79.9  kcal  mol  1  for 
the  theoretical  calculations  and  (ECou)  =  83  kcal  mol  1  for  the  experiment).  Isotopic  substitution 
was  used  in  the  experiment  to  distinguish  between  the  reactive  channel,  160  +  C180  — >  l6OC  + 
l80,  Eq.  (2),  and  the  inelastic  channel,  160  +  C180  — ►  160  +  C180,  Eq.  (1).  We  investigated  the 
dynamics  of  the  reactive  and  inelastic  pathways,  including  the  disposal  of  energy  in  the  products 
and  the  angular  distributions  of  the  scattered  products.  Using  the  surfaces  of  Ref.  7,  new 
quasiclassical  trajectory  (QCT)  calculations  were  performed  which  extend  previous  results  to 
facilitate  direct  comparisons  with  the  present  measurements  and  to  more  fully  explore  the 
dynamics.  The  results  of  the  QCT  calculations  are  in  almost  quantitative  agreement  with  the 
present  experimental  data. 
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II.  THEORETICAL  DETAILS 


A.  Methods.  The  theoretical  methods  are  based  on  the  potential  surfaces  and  quasi- 
classical  trajectory  (QCT)  procedures  of  Braunstein  and  Duff.7  We  use  the  same  surfaces  of 
Braunstein  and  Duff,  and  we  extend  the  QCT  calculations  to  generate  new  cross  section  results 
that  can  be  compared  directly  with  the  experimental  data  discussed  in  this  paper  and  at  several 
collision  energies  well  below  and  above  those  used  in  the  experiment.  The  present  work 
provides  new  and  extensive  results  that  separate  the  reactive  and  inelastic  cross  sections.  The 
present  results  also  use  at  least  an  order  of  magnitude  more  trajectories  per  collision  energy  than 
the  earlier  work  for  improved  statistical  convergence. 

Potential  energy  surfaces.  We  briefly  review  the  key  features  of  the  potential  energy 
surfaces  used.  Complete  details  can  be  found  elsewhere.7  Figure  1  shows  the  basic  features  of 
the  relevant  potentials.  Collisions  of  0(  P )  with  CO(  X  )  correlate  to  three  electronic  states  of 
CO2,  all  triplets:  a  single  A'  state  and  two  A"  states.  The  five  nearby  singlet  states,  which  can 
only  participate  through  weak  spin-orbit  interactions,  lead  to  Of1/))  and  CO(*X+),  about  46  kcal 
moT1  above  the  Of3/5)  +  COf  *X  )  threshold.  One  of  these  five  singlet  states  leads  to  the  *Xg+ 
ground  state  of  CO2  about  127  kcal  moT1  below  the  Of3/5)  +  COf  *X  )  threshold.  The  remainder 
are  relatively  high  in  energy  and  give  rise  to  some  of  the  complex  features  in  the  CO2  UV 
photoabsorption  spectrum  starting  about  23  kcal  mol  1  above  the  Of3/5)  +  COf'X  )  threshold.11 
As  in  our  earlier  work,7  we  ignore  spin-orbit  interactions  that  could  mix  the  singlet  and  triplet 
states. 

The  potentials  used  are  global  fits  of  -320  separate  ab  initio  calculations  of  the  three 
lowest  triplet  states  of  CO2.  The  ab  initio  calculations  were  performed  at  the  (12  electron,  10 
orbital)  CASSCF-MP2  level,  within  a  modest  631+Gfd)  basis  set,  with  the  electronic  structure 
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code,  GAMESS.12  The  lowest  two  states,  3 A'  and  1  3A",  have  well  depths  of  21.7  kcal  mol  1  and 
25.1  kcal  mol  ',  respectively.  These  minima  occur  near  the  geometries  shown  in  Figure  1.  For 
the  3 A'  state  the  minimum  occurs  at  ri=r2=1.26  A,  9  =  118,  and  for  the  1  3 A"  state  the  minimum  occurs 
at  ri=r2=  1.27  A,  9  =  127°  The  2  3 A"  state  is  mostly  repulsive  and  has  a  saddle  point  with  C2v 
geometry.  The  3 A'  and  1  3 A "  states  have  small  barriers  of  4.6  kcal  mol  1  and  6.9  kcal  mol-1, 
respectively,  which  occur  at  Cs  geometries  of  their  transition  states  near  to  those  shown  in  Figure 
1.  These  transition  states  are  in  bent  Cs  geometries,  where  for  the  3 A'  state,  r\  =  1.16  A,  r2  =  2.9 
A,  9  =  1 12°,  and  for  the  1  3 A"  state,  r\  =  1.16  A,  r2  =  1.85  A,  9  =  122°.  The  ab  initio  points  of 

13 

the  three  electronic  states  were  globally  fit  with  the  method  of  Aguado  and  Paniagua. 

Dynamics  and  cross  sections.  The  three  lowest  triplet  potential  energy  surfaces  were 
used  to  generate  cross  sections  at  several  fixed  collision  energies.  Each  potential  surface  was 
treated  independently,  and  total  cross  sections  were  obtained  by  weighting  the  contribution  of 
each  electronic  state  by  1/3:  <j=  l/3[a(3A')  +  a(l  3 A")  +  a(2  3 A")].  We  used  standard  Monte 
Carlo  methods14  to  generate  quasiclassical  trajectory  (QCT)  cross  sections,  separating 
contributions  into  reactive  (r)  collisions,  where  the  incoming  O  atom  is  exchanged  with  the  O 
atom  belonging  to  the  target  CO,  and  inelastic  (i)  collisions,  where  the  incoming  O  atom  is  not 
exchanged  with  the  O  atom  belonging  to  the  target  CO.  For  a  trajectory  to  contribute  to  the 
inelastic  cross  section,  the  final  CO(v',j')  state  must  be  different  than  the  initial  CO(v,j)  state. 
We  note  that  all  three  electronic  states  will  have  significant  contributions  to  the  inelastic  cross 
sections.  The  2  A"  state  will  in  general  have  only  minor  contributions  to  the  reactive  cross 
sections  because  of  its  repulsive  nature.  The  lowest  two  electronic  states,  3A'  and  1  3 A",  will 
have  major  contributions  to  the  reactive  cross  sections  because  of  their  low  barriers  to  reaction. 
All  calculations  were  performed  with  12C  and  160  atom  masses  for  the  target  CO  and  l60  for  the 
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incoming  O-atom.  The  CO  mass  used  in  the  calculations  is  therefore  slightly  different  from  the 
experimental  measurements,  which  used  160  +  12C180.  However,  at  relative  collision  velocities 
of  ~8  km  s  ,  we  believe  that  differences  in  results  arising  from  the  mass  difference  between 
12C160  and  12C180  are  much  smaller  than  other  approximations  made  in  the  present  calculations. 
All  calculations  were  done  with  the  target  CO  in  the  ground  vibrational  state  and  a  rotational 
temperature  of  300  K. 

We  performed  two  sets  of  QCT  calculations.  In  one  set  we  focused  on  obtaining  results 
that  could  be  directly  compared  to  the  present  measurements.  We  set  the  collision  velocity  to  8 
km  s_1  (77.9  kcal  mol-1).  The  maximum  impact  parameter  was  set  to  10  a.u.  (5.29  A),  and  1><106 
trajectories  were  run  for  each  of  the  three  electronic  states.  The  value  of  the  maximum  impact 
parameter  was  made  relatively  large  in  order  capture  the  low-angle,  low-energy  (rotational) 
scattering  in  the  inelastic  channel,  which  is  important  when  comparing  to  the  present 
measurements.  We  calculated  differential  angular  cross  sections, 

da  _  NAe  1 

dQ,  Ntot  2n(sm0)k0  ’ 

where  NAe  is  the  number  of  trajectories  within  the  angular  bin,  AO,  and  Nm  is  the  total  number  of 
trajectories.  To  compare  with  the  fine  angular  resolution  of  the  present  measurements,  we  set 
A0=  2°.  As  discussed  in  detail  later,  for  the  reactive  channel,  the  c.m.  scattering  angle  0  is  taken 
to  be  the  angle  between  the  c.m.  velocity  vector  of  the  reagent  O  atom  and  the  c.m.  velocity 
vector  of  the  product  CO.  For  the  inelastic  channel,  0  is  the  angle  between  the  c.m.  velocity 
vector  of  the  reagent  O  atom  and  the  c.m.  velocity  vector  of  the  product  O  atom.  We  also 
calculated  c.m.  differential  translational  energy  cross  sections, 
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da 


N 


AEtrans 


1 


dE„ 


AC,  AE„ 


where  j\[A,;t,ans  js  the  number  of  trajectories  with  a  final  product  translational  energy,  Etrans,  within 
an  energy  bin,  AEtram-  We  used  an  energy  bin  width  of  5.8  kcal  mol  L,  which  approximates  the 
experimental  translational  energy  resolution. 

In  a  second  set  of  QCT  calculations,  we  focused  on  the  collision  energy  dependence  of 
the  rovibrationally  resolved  cross  sections.  We  performed  QCT  calculations  at  O  +  CO  collision 
velocities  of  3,  4,  5,  6,  7,  8,  9,  and  10  km  s-1,  corresponding  to  c.m.  collision  energies  of  11.0, 
19.5,  30.4,  43.8,  59.6,  77.9,  99.6,  and  121.7  kcal  mol  respectively,  with  the  maximum  impact 
parameter  set  to  5  a.u.  (2.65  A).  The  maximum  impact  parameter  of  5  a.u.  used  in  these 
calculations  will  not  be  adequate  to  capture  low-angle  inelastic  scattering  resulting  in  CO(v=0) 
products  (pure  rotational  transitions).  However,  these  calculations  should  be  sufficient  for  all 
reactive  collisions  and  inelastic  collisions  leading  to  CO  vibrational  excitation.  In  this  set  of 
calculations,  2.5  x  105  trajectories  were  run  for  each  electronic  state  at  each  collision  energy. 
This  number  of  trajectories  is  about  a  factor  of  10  larger  than  used  in  our  earlier  study.7  In  all 
calculations,  standard  histogram  binning  was  used. 

B.  Theoretical  Results,  Results  for  8  km  s~‘  collision  velocity.  Figure  2  shows  the 
calculated  results  for  the  differential  translational  energy  cross  sections  for  reactive,  Eq.  (2),  and 
inelastic,  Eq.  (1),  channels  of  O  +  CO  collisions  at  a  relative  velocity  of  8  km  s_1  (77.9  kcal  mol” 
').  The  reactive  and  inelastic  distributions  are  shown  with  (solid  lines)  and  without  (dashed 
lines)  the  contribution  of  product  CO(v=0),  to  separate  the  dynamics  of  pure  rotational 
excitation,  which  dominates  the  inelastic  channel  at  high  translational  energies,  from  vibrational 
excitation.  The  energy  bin  width  is  equal  to  5.8  kcal  mol-1  starting  at  a  center  bin  of  2.9  kcal 
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mol-1,  so  that  these  theoretical  results  are  directly  comparable  to  the  present  measured 
distributions,  which  have  a  similar  resolution.  Distributions  obtained  after  quantum-binning  the 
product  CO(v,j)  states,  discussed  further  below,  should  reveal  detailed  structure  corresponding  to 
rovibrationally  resolved  CO.  We  note  that  the  highest  energy  non-zero  cross  sections  extend  to 
the  translational  energy  bin  centered  on  83.7  kcal  mol-1,  which  is  above  the  collision  energy. 
Because  the  reagent  CO  is  at  300  K,  some  initial  rotational  energy  ends  up  in  product  translation 
for  a  small  fraction  of  the  collisions,  and  there  is  some  leakage  into  this  high  energy  bin. 

In  Fig.  2,  a  large  product  translational  energy,  Ctrans,  corresponds  to  a  low  product  internal 
energy.  The  peak  and  general  behavior  at  high  F'trans  for  the  inelastic  channel  corresponds  to 
mostly  pure  rotational  excitation,  and  it  is  typical  of  non-reactive  collisions  that  occur  on  a 
repulsive  potential.  At  lower  /7,-ans  (higher  CO  product  internal  energy),  particularly  below  £n-ans 
=  40  kcal  mol  ,  the  reactive  and  inelastic  channels  have  the  same  broad  and  flat  shape  and 
nearly  the  same  magnitude.  The  appreciable  value  of  the  differential  cross  section  persists  to 
near  the  zero  translational  energy  cut-off,  corresponding  to  all  the  available  collision  energy 
going  into  CO  internal  energy.  Also  shown  are  the  reactive  and  inelastic  differential 
translational  energy  cross  sections  with  the  contribution  from  the  product  CO(v=0)  removed — 
i.e.,  all  trajectories  where  the  product  CO  is  in  the  vibrational  ground  state  do  not  contribute. 
The  reactive  and  inelastic  channels  with  the  CO(v=0)  contribution  removed  are  nearly  identical 
for  all  energies  and  track  the  full  reactive  channel  (including  the  CO(v=0) 
contribution  (differential  cross  sections  below  about  Claris  ~  70  kcal  mol  ,  which  is  the 
approximate  threshold  for  creating  vibrationally  excited  CO.  This  result  implies  that  the 
underlying  dynamics  for  collision  products  where  CO  is  vibrationally  excited  are  similar  for  both 
the  reactive  and  inelastic  channels. 
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Figure  3  shows  the  calculated  differential  angular  cross  section  for  O  +  CO  collisions  at  8 
km  s-1  collision  velocity  for  the  reactive,  Eq.  (2),  and  inelastic  channels,  Eq.  (1).  Results  are 
shown  both  with  (solid  lines)  and  without  (dashed  lines)  the  product  CO(v=0)  contribution,  in 
order  to  separate  pure  rotational  excitation  from  rovibrational  excitation.  Use  of  a  relatively 
small  angular  bin  width  of  2  degrees  reveals  small  rapid  oscillations  in  these  curves  for  smaller 
values  of  the  cross  sections,  which  is  an  indication  of  the  statistical  uncertainty  of  these  results. 
In  Fig.  3,  the  large  values  at  low  angles  for  the  inelastic  channel  are  typical  of  low-angle,  low- 
energy-transfer  collisions  on  repulsive  surfaces.  The  peak  at  low  angle  for  the  reactive  channel 
suggests  a  “stripping”  mechanism,  where  the  incoming  O-atom  strips  off  a  carbon  atom  from  the 
target  CO,  and  the  product  CO  is  scattered  mostly  forward — i.e.,  in  the  same  direction  as  the 
reagent  O-atom  velocity  vector.  Whether  it  reacts  or  is  scattered  inelastically,  the  reagent  O 
atom  is  scattered  mostly  in  the  same  direction  as  its  initial  velocity  vector.  The  transition  state 
for  reaction  is  bent,  near  120°,  which  may  help  explain  how  the  CO  reaction  product  scatters  in 
the  same  direction  as  the  initial  direction  of  the  incoming  O  atom.  The  differential  angular  cross 
section  for  the  reactive  channel  is  broad  as  a  function  of  scattering  angle,  with  a  small  backward 
peak  near  180°.  Also  shown  are  the  differential  angular  cross  sections  without  the  CO(v=0) 
contributions.  The  reactive  and  inelastic  cross  sections  without  the  CO(v=0)  contributions  are 
again  very  similar  and  closely  approximate  the  full  reactive  differential  cross  section.  It  may  be 
that  for  vibrational  excitation,  reactive  and  inelastic  trajectories  follow  similar  paths,  either  on 
the  reagent  side  of  the  transition  state  barrier  for  inelastic  collisions  or  product  side  of  the 
transition  state  barrier  for  reactive  collisions  (see  Fig.  12  of  Braunstein  and  Duff7). 

Figures  4a-c  show  the  reactive  (Fig.  4a)  and  inelastic,  with  (Fig.  4b)  and  without  (Fig.  4c) 
the  product  CO(v=0)  contribution,  rovibrationally  resolved  cross  sections  versus  the  product 
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translational  energy.  The  values  of  Etmns  in  these  plots  were  obtained  by  Etmns  =  77.9  kcal  mol  1 
-  £intemai(CO(v,j)).  The  numbers  in  Figs.  4a-c  (0-7)  correspond  to  CO(v)  level  thresholds.  (The 
<j(v,j=0)  cross  sections  dip  down  to  very  low  values  and  so  conveniently  indicate  vibrational 
level  onsets.)  For  the  inelastic  cross  sections,  Figs.  4b  and  4c,  much  of  the  v=0  cross  section  is 
above  the  scale  of  the  figures.  This  expression  for  Tuans  neglects  a  small  energy  contribution 
from  the  initial  reagent  CO  internal  energy.  We  note  also  that  the  cross  section  is  obtained  after 
binning  the  classical  trajectories  into  quantum  states.  Therefore,  these  finely  resolved  cross 
sections  cannot  be  compared  quantitatively  to  the  differential  translational  energy  cross  sections 
shown  in  Fig.  2.  However,  they  are  nearly  the  same  quantity,  and  they  reveal  rich  and  detailed 
state  populations.  As  we  shall  show,  product  CO  vibrational  excitation  for  v>10  and  j>100  is 
seen  in  the  calculations.  Many  of  these  CO(v,j)  states  are  evident  in  Fig.  4,  but  many  of  the 
higher  (v,j)  states  overlap  the  thresholds  of  higher  vibrational  manifolds,  making  these  states 
difficult  to  distinguish  from  one  another.  This  rich  structure  resulting  from  the  highly  excited 
state  populations  shown  in  Fig.  4  underlie  the  results  of  Fig.  2  and  could  possibly  be  revealed  in 
very  highly  resolved  energy  measurements.  As  seen  in  Fig.  2,  the  overall  magnitude  and  shape 
of  the  inelastic  rovibrationally  resolved  cross  sections  without  the  product  (CO(v=0)) 
contribution,  Fig.  4c,  closely  resembles  the  full  reactive  rovibrational  distribution,  Fig.  4a.  The 
largest  differences  occur  at  extremely  low  Tlrans  (high  CO  internal  energy)  where  the  reactive 
cross  sections  remain  fairly  large  up  until  the  Etmns  =  0  cut-off. 

Figure  5  shows  the  reactive  (Fig.  5a)  and  inelastic  (Fig.  5b)  rovibrationally  resolved  cross 
section  a(CO(v,j)),  for  O  +  CO  collisions  at  8  km  s  for  v=0-3.  We  have  multiplied  the  a(0,j) 

_3 

inelastic  cross  sections  by  1x10  to  put  these  data  on  a  common  scale.  The  rapid  oscillations  in 
the  other  curves  are  an  indication  of  the  statistical  uncertainty  in  the  Monte  Carlo  calculations. 
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The  reactive  cross  sections  for  each  vibrational  level  have  similar  broad  shapes,  with  single 
maxima,  that  extend  to  the  maximum  available  energy,  except  for  a(0,j)  which  has  two  maxima, 
one  at  high  j  and  one  at  low  j.  For  cr(v>0,j)  taken  together,  the  reactive  rotational  distributions 
can  be  fit  well  with  a  Boltzmann  rotational  distribution  of  ~8,000  K.  (The  vibrational 
temperature  is  also  near  8,000  K.)  The  rovibrationally  resolved  cross  sections  persist  to  very 
high  rotational  quantum  numbers,  up  to  j~l  1 0  (corresponding  to  ~65  kcal  mol  ),  which  is  an 
appreciable  fraction  of  the  collision  energy  of  77.9  kcal  mol  The  inelastic  rovibrationally 
resolved  cross  sections  for  product  CO(v>0)  have  a  similar  broad  shape  as  the  reactive  cross 
sections,  extending  to  j~100. 

Results  at  collision  velocities  from  3  km  s~‘  through  10  km  s  1 .  To  gain  more  insight  into 
the  dynamics,  and  in  anticipation  of  possible  future  experiments,  Fig.  6  shows  the  vibrationally 
resolved  cross  sections,  for  O  +  CO  reactive  and  inelastic  collisions  at  3,  4,  5,  6,  7,  8,  9,  and  10 
km  s"1  collision  velocities.  The  inelastic  a(v=0)  cross  section  is  above  the  scale  of  the  figure. 
Except  for  the  highest  vibrational  populations  near  the  available  energy  limit  for  a  particular 
collision  velocity,  the  reactive  vibrational  distributions  are  similar  to  Boltzmann  distributions. 
The  effective  temperatures  of  these  distributions  increase  with  approximately  the  square  of  the 
collision  velocity.  The  inelastic  vibrationally  resolved  cross  sections  closely  resemble  the 
reactive  vibrationally  resolved  cross  sections  except  for  a(v=0).  For  the  present  energy 
dependent  results,  which  focus  on  CO  vibrational  excitation  and  where  the  maximum  impact 
parameter  is  5  a.u.,  the  ct(v=0)  inelastic  results  will  not  include  significant  contributions  from 
large  impact  parameter,  small  A J  transitions.  For  vibrationally  excited  products,  the  inelastic  CO 
distributions  are  slightly  hotter  than  the  reactive  distributions. 
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Figure  7  examines  the  average  rotational  and  vibrational  energies  of  product  CO,  where 
we  show  the  rotational  and  vibrational  energy  fractions,  /J  and  /v  respectively,  for  O  +  CO 
reactive  and  inelastic  collisions  at  3,  4,  5,  6,  7,  8,  9,  and  10  km  s_1  collision  velocities.  The 
rotational  and  vibrational  energy  fractions  are  defined: 


fJ  = 


J  collision 


r  = 


\  ^ collision 


J  collision 


jaj 


J  collision 


7-Z £ . 


respectively,  where  Ej  is  the  rotational  energy,  EY  is  the  vibrational  energy,  Ecomsion  is  the 

collision  energy,  <jtot  is  the  total  cross  section,  <j(j)  =  I  0"(v,y),  and  ofv)  ^^rr(v,y).  The 

v  j 

vibrational  energy  fraction  increases  with  collision  energy  for  the  reactive  and  inelastic  channels. 
The  reactive  energy  fraction  shows  a  rapid  increase  from  threshold,  which  may  reflect  the  small 
energy  barrier  to  reaction.  The  large  o(v=0)  inelastic  cross  section  makes  the  overall  inelastic 
vibrational  energy  fraction,  which  is  an  average  of  the  energy  of  vibrational  states,  smaller  than 
the  reactive  energy  fraction  The  reactive  rotational  energy  fraction  is  generally  flat  while  the 
inelastic  rotational  energy  fraction  decreases  with  collision  energy.  For  collision  velocities 
above  8  km  s  ,  the  distribution  of  internal  energy  in  the  elastic  channel  continues  to  change,  but 
it  remains  relatively  constant  in  the  reactive  channel. 

Figure  8  shows  two  different  excitation  functions  for  O  +  CO  collisions:  a  photon 
excitation  function,  <jp,  and  a  vibrational  excitation  function,  <je,  both  as  a  function  of  collision 
energy.  These  excitation  functions  are  separated  into  reactive  and  inelastic  components.  We 
define  the  photon  excitation  function  as  ap  =^/cr(v  =  /) .  The  photon  excitation  function  is 

i=i 
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particularly  convenient  for  comparisons  to  laboratory  measurements  of  photon  production  due  to 
O  +  CO  collisions  which  integrate  signal  over  the  entire  fundamental  CO  band  from  4-7  microns, 
such  as  the  measurements  of  Upschulte  and  Caledonia.5  The  photon  excitation  function  also 
provides  a  way  to  compare  directly  reactive  and  inelastic  excitation  functions  by  not  including 
CO  products  formed  in  the  vibrational  ground  state,  which  for  the  inelastic  channel  would 
otherwise  mask  the  behavior  of  the  vibrationally  excited  CO  cross  section.  The  photon 
excitation  functions  for  the  reactive  and  inelastic  channels  are  nearly  the  same,  with  the  inelastic 
channel  slightly  larger.  This  behavior  reflects  the  similarity  of  the  underlying  vibrationally 
resolved  cross  sections  for  excited  CO  products.  We  define  the  vibrational  excitation  function 
as  cje  =  y  a(  v  =  i) .  The  vibrational  excitation  functions  are  also  very  similar. 

;=i 


III.  EXPERIMENTAL  DETAILS 

A.  Methods.  The  experiments  were  performed  with  the  use  of  a  crossed-molecular 
beams  apparatus  (see  Fig.  9)  equipped  with  a  fast-atom  beam  source.8'15’18  General  details  about 
the  experimental  apparatus  and  details  of  the  analysis  can  be  found  in  an  earlier  paper.  A  pulsed 
beam  of  oxygen  atoms  (160)  was  crossed  at  right  angles  with  a  pulsed,  supersonic  beam  of 
C  O  gas.  Products  that  scattered  from  the  interaction  region  were  detected  with  a  rotatable 
mass  spectrometer  detector  that  measured  number-density  distributions  as  a  function  of  arrival 
time  at  the  electron-bombardment  ionizer,  N(t).  These  number-density  distributions  are 
commonly  referred  to  as  time-of-flight  (TOF)  distributions.  Laboratory  angular  distributions, 
N(0),  where  0  is  the  laboratory  angle  of  the  scattered  products  with  respect  to  the  O-atom  beam, 
are  determined  by  integrating  the  TOF  distributions  over  time  at  a  fixed  angle  for  a  series  of 
laboratory  angles,  0.  These  data  are  collected  in  the  laboratory  reference  frame,  so  a  forward 
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convolution  method19'20  is  employed  to  derive  center-of-mass  (c.m.)  quantities  from  the 
laboratory  TOF  and  angular  distributions. 

The  hyperthennal  beam  used  in  the  experiment  had  an  average  translational  energy  of 
115  kcal  mol  1  (Fig.  10),  and  the  mole  fraction  of  atomic  oxygen  in  the  beam  was  approximately 
98  percent.  The  hyperthennal  160-atom  beam  was  crossed  at  a  90°  intersection  angle  by  a  pulsed 

beam  of  isotopically-labeled  carbon  monoxide  molecules  (C  O).  The  C  O  (Isotec,  Sigma 

18  16 

Aldrich)  was  98.7%  C  O,  with  the  remainder  (1.3%)  being  C  O.  A  supersonic  expansion  of 

18  18 
C  O  was  created  with  a  piezoelectric  pulsed  valve  from  a  stagnation  pressure  of  10  psig  C  O. 

The  C  O  beam  passed  through  a  5 -nun-diameter  skimmer  and  a  3 -nun-diameter  aperture  before 

crossing  the  O-atom  beam.  The  distance  between  the  nozzle  and  skimmer  was  9.2  cm,  and  the 

distance  from  the  skimmer  to  the  3-mm-diameter  aperture  was  2.8  cm.  From  the  aperture,  the 

beam  traveled  1.5  cm  to  the  crossing  point  of  the  two  beams.  The  C  O  beam  velocity  was 

estimated  to  be  800  ±  80  m  s'1.21  The  velocity  of  the  C180  beam  was  an  order  of  magnitude 

lower  than  the  O-atom  beam,  and  the  velocity  width  of  the  C  O  beam  was  not  considered  in  the 

analysis  of  the  experimental  data.  The  crossing  region  of  the  two  beams  was  located  99  cm 

from  the  apex  of  the  conical  nozzle  of  the  O-atom  source  and  13  cm  from  the  orifice  of  the  C  O 

pulsed  valve  source. 

The  collision-energy  distribution  in  the  c.m.  reference  frame  was  derived  from  the  O- 
atom  beam-velocity  distribution  (related  to  Fig.  10  by  Etr.u^  =  \  mv2 )  and  the  nominal  velocity  of 

the  carbon  monoxide  beam  using  the  equation,  Ecoll  =  \juv2rel ,  where  pi  is  the  reduced  mass  and 

Vrei  is  the  relative  velocity  of  the  O  atoms  and  the  C  O  molecules.  The  average  collision  energy 
in  this  experiment  was  83  kcal  mol  \  and  the  width  of  the  collision-energy  distribution  (FWHM) 
was  ~30  kcal  mol1,  ranging  from  70  to  100  kcal  mol  \  A  Newton  diagram  for  collisions  of  O 
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atoms  and  C180  with  Eco\\  =  83  kcal  mol  1  is  shown  in  Fig.  11.  This  diagram  shows  maximum 
recoil  velocities  for  inelastically-scattered  C  O  (solid  blue  circle)  and  O  (dashed  blue  circle). 
It  also  shows  maximum  recoil  velocities  for  reactively-scattered  OC  (solid  red  circle)  and  O 
(dashed  red  circles).  The  shaded  region  shown  on  the  Newton  diagram  indicates  the  range  of 
laboratory  angles  used  for  the  experiment:  6°  to  54°  with  respect  to  the  nominal  velocity  vector 
of  the  hyperthennal  O-atom  beam. 

The  experiments  focused  on  the  reactively-scattered  16OC  molecules  (m/z  =  28), 
inelastically-scattered  C  O  molecules  ( m/z  =  30),  and  inelastically-scattered  O  atoms  ( m/z  =  16). 
Signals  from  m/z  =  32  (02+)  were  used  to  correct  the  m/z  =16  data  for  the  contribution  of 
inelastically  scattered  O2  cracking  to  0+  on  the  ionizer.  Low  signals  combined  with  high 
background  levels  of  water  in  the  detector  at  m/z  =18  made  it  impossible  to  detect  the  reactive 
l80  product.  TOF  distributions  were  collected  for  2000  beam  pulses  for  m/z  =  28  (16OC+)  and  30 
(C180+),  and  200  beam  pulses  for  16  (0+),  and  32  (C>2+)  at  each  laboratory  angle.  The  laboratory 
angle  was  adjusted  in  2  degree  increments  until  the  entire  angular  range  was  covered.  Then  the 
increment  direction  was  reversed  and  the  cycle  was  repeated  until  a  total  of  four  TOF 
distributions  had  been  collected  for  each  product  at  each  detector  angle.  Because  this  process 
lasted  more  than  35  hours,  it  was  important  to  sum  the  TOF  distributions  collected  in  this  manner 
to  account  for  minor  long-term  drifts  in  the  experimental  parameters.  At  a  representative 
laboratory  angle  of  10°,  the  integrated  count  rates  were  8.5><106  counts  s_1  for  m/z  =  30  (C180+), 
1.2*108  counts  s_1  for  m/z  =  16  (0+),  and  2.5 xlO7  counts  s^1  for  m/z  =  28  (16OC+). 

B.  Experimental  Results.  Inelastic  scattering.  Figure  12  shows  TOF  distributions  for 
m/z  =  16  (0+)  and  30  (ClsO+)  collected  at  five  representative  laboratory  detector  angles.  Figure 
13  shows  laboratory  angular  distributions  for  inelastically  scattered  O  and  C  O.  The  curves  in 
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Figs.  12  and  13  are  the  forward-convolution  fits  to  the  laboratory  angular  distributions  for  the  O 

product  (red)  and  the  C180  product  (blue),  based  on  the  can.  translational  energy  and  angular 

distributions  in  Fig.  14.  The  distributions  shown  in  Fig.  14  are  taken  from  the  theoretical  results 

in  Figs.  2  and  3  but  presented  on  a  linear  scale.  These  calculated  c.m.  distributions  accurately 

predict  the  laboratory  distributions,  as  seen  in  Figs.  12  and  13.  The  O  product  is  forward 
18. 

scattered  and  the  C  O  is  backward  scattered  with  respect  to  the  initial  direction  of  the  reagent  O 

atoms  in  the  c.m.  frame.  The  laboratory  detection  angles  are  primarily  sensitive  to  C  O 

molecules  that  scatter  into  the  backward  hemisphere  in  the  c.m.  reference  frame  with  respect  to 

the  direction  of  the  incident  O  atoms  (see  Fig.  1 1),  so  only  the  angular  range  in  the  c.m.  angular 

distributions  that  was  accessible  in  experiment  was  used  in  the  forward  convolution.  The 

experiments  were  sensitive  to  C  O  and  O  scattered  from  0°  to  120°  and  from  0°  to  80°, 

respectively,  in  the  c.m.  reference  frame.  The  translational  energy  distribution  reveals  that,  on 

average,  about  85%  of  the  available  energy  (83  kcal  mol  ')  is  released  in  translation  in  the 

inelastic  channel,  although  there  is  a  tail  to  low  translational  energies  suggesting  that  some 

collisions  result  in  large  energy  transfers  to  internal  degrees  of  freedom  of  C  O.  The  angular 

18. 

distributions  of  scattered  O  and  C  O  illustrate  the  propensity  for  forward  scattering  which 

results  from  high-impact-parameter  collisions. 

Although  the  hyperthennal  beam  in  our  experiment  had  only  two  percent  CF,  the  inelastic 

cross  section  for  O2  +  C  O  is  almost  twice  that  for  O  +  C  O,  so  the  contribution  to  the  inelastic 

scattering  signal  detected  at  m/z  =  30  (C  O  )  was  not  beneath  our  detection  limits.  A  second 

18  18 

inelastic  scattering  process,  O2  +  C  O  — >  C  O  +  O2  was  identified  as  a  small,  fast  shoulder  in 
the  TOF  distributions  collected  at  m/z  =  30  (C  O  ).  A  functional  fonn  for  the  fast  component 
was  determined,  and  this  component  was  subtracted  from  all  the  TOF  distributions  collected  at 
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18  + 

m/z  =  30  (C  O  ),  including  those  shown  in  Fig.  12.  Only  the  slower  (main)  peak,  corresponding 
to  the  C  O  that  scattered  inelastically  from  O,  was  used  to  gather  dynamical  infonnation  about 
the  inelastic  channel. 

Reactive  scattering.  TOF  distributions  collected  at  m/z  =  28  (16OC+)  at  five 
representative  laboratory  angles  are  shown  in  Fig.  15,  and  the  corresponding  laboratory  angular 
distribution  is  shown  in  Fig.  16.  Center-of-mass  angular  and  translational  energy  distributions 
for  the  0(3P )  +  C180  — >  16OC  +  180  reaction  were  obtained  from  the  QCT  calculations  and  are 
shown  as  dashed  blue  curves  in  Fig.  17.  These  c.m.  distributions  were  used  in  the  forward- 
convolution  procedure  to  predict  laboratory  TOF  and  angular  distributions  shown  as  dashed  blue 
curves  in  Figs.  15  and  16.  The  TOF  distributions,  predicted  from  the  theoretical  c.m. 
translational  energy  and  angular  distributions,  match  the  experimental  distributions  very  well, 
but  the  predicted  laboratory  angular  distribution  does  not  quite  match  the  observed  flux 
difference  between  forward  and  sideways  scattering  in  the  laboratory  angular  distribution.  The 
experimentally  derived  c.m.  translational  energy  and  angular  distributions  are  shown  as  red 
curves  in  Fig.  17.  These  experimental  curves  were  obtained  by  optimizing  a  translational  energy 
distribution,  based  on  an  RRK  form”  ,  and  an  angular  distribution  of  arbitrary  point  form.  The 
corresponding  laboratory  TOF  and  angular  distributions  are  shown  as  the  solid  red  curves  in 
Figs.  15  and  16.  The  theoretically  calculated  translational  energy  distribution  exhibited  some 
structure,  but  the  experimental  resolution  did  not  justify  adding  structure  to  the  experimentally 
derived  translational  energy  distribution.  Thus,  within  experimental  error,  both  the  theoretical 
and  experimental  translational  energy  distributions  may  be  the  same.  The  difference  between  the 
theoretical  and  experimental  c.m.  angular  distributions  is  believed  to  be  just  outside  experimental 
error.  The  theoretical  and  experimental  results  both  indicate  that  about  60%  of  the  available 
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energy  is  released  in  translation  and  that  the  16OC  product  is  predominantly  forward  scattered. 
The  theoretical  calculations  predict  a  small  increase  in  the  amount  of  16OC  scattering  at  the  most 
backward  angles,  near  180°,  but  that  angular  range  is  beyond  the  range  to  which  the  experiment 
is  sensitive.  A  velocity-flux  map,  shown  in  Fig.  18,  illustrates  the  differential  scattering  cross 
section  of  the  reactive  16OC  product,  where  16OC  tends  to  be  forward  scattered,  with  a  broad 
range  of  product  velocities.  In  the  angular  range  probed  by  the  experiment,  the  theoretical  and 
experimental  results  for  this  channel  are  in  almost  quantitative  agreement,  with  only  a  slight 
discrepancy  in  the  angular  distribution. 

IV.  DISCUSSION 

The  experimentally  observed  inelastic  scattering  data  is  predicted  quantitatively  by  the 
theoretical  c.m.  translational  energy  and  angular  distributions.  As  expected  for  inelastic 
collisions,  O  atoms  scatter  from  CO  with  little  change  in  direction  and  velocity.  Thus,  the  O 
atoms  are  predominantly  forward  scattered  and  the  CO  molecules  are  predominately  backward 
scattered  in  the  c.m.  reference  frame.  And  the  c.m.  translational  energy  distribution  has  a 
maximum  near  the  c.m.  collision  energy.  Subtracting  the  translational  energy  distribution  from 
the  collision  energy  gives  the  internal  energy  distribution  of  scattered  CO.  On  average,  -15 
percent  of  the  collision  energy  is  transferred  into  internal  energy.  However,  the  low-energy  tail 
in  the  translational  energy  distribution  indicates  that  a  large  fraction  (-84  percent)  of  single 

3  1 

o  CP)  +  CO  collisions  can  transfer  more  than  —6  kcal  mol  into  internal  energy  in  the  CO, 
corresponding  to  greater  than  one  quanta  of  vibration  (v=l).  A  smaller  fraction  (~0.5  percent)  of 
scattered  CO  molecules  may  have  more  than  -63  kcal  mol  1  of  internal  energy,  corresponding  to 
CO  (v=ll).  Although  the  majority  of  inelastic  collisions  lead  to  CO  that  is  scattered  in  the 
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backward  direction  with  respect  to  the  initial  direction  of  the  reagent  O  atoms,  it  should  be  noted 
that  C180  detected  in  the  experiment  can  only  come  from  collisions  of  O  atoms  with  C180 
molecules  that  lead  to  C  O  scattering  in  the  forward  direction  with  respect  to  the  reagent  O 
atoms  (refer  to  Newton  diagram  in  Fig.  11).  The  forward  scattered  CO  could  be  the  result  of 
low-impact-parameter  collisions.  Large  energy  transfers  have  been  observed  when  hyperthermal 
Ar  atoms  scatter  from  C2H6  or  C2F6  at  low  impact  parameters;  ’  however,  in  analyzing  the 
calculated  O  +  CO  trajectories,  it  was  difficult  to  conclude  that  low-impact-parameter  collisions 
were  necessarily  correlated  with  large  energy  transfers.  Instead,  large  energy  transfers  in 
inelastic  collisions  seem  to  be  correlated  with  failed  reactive  collisions  (see  below). 

An  analysis  of  the  experimental  results  suggests  that  roughly  ten  percent  of  the  collisions 
of  0(  P)  with  CO  at  (Econ)  =  83  kcal  mol  lead  to  an  O-atom  exchange  reaction,  detected  by  the 
observation  of  scattered  16OC.  Assuming  a  typical  hard-sphere  cross  section  of  ~10  15  cm2,  this 
experimental  result  would  suggest  a  total  reactive  cross  section,  ^ <jr(y  =  i ),  of  -TO  16  cm2. 

i= 0 

This  value  is  in  rough  agreement  with  the  present  QCT  calculated  total  reactive  cross  section  of 

1.98x10'  cm-.  This  rough  agreement  would  seem  to  support  the  conclusion  that  the  vibrational 

excitation  cross  section  magnitudes  reported  earlier  by  Upschulte  and  Caledonia5  and  Green  et 

a/.6  are  too  low  by  about  an  order  of  magnitude.  Within  experimental  uncertainty,  the 

theoretically  calculated  and  experimentally  observed  c.m.  translational  energy  distributions  for 

reactive  collisions  (see  Fig.  17)  are  identical,  suggesting  a  broad  range  of  product  translational 

energies  with  an  average  of  60  percent  of  the  available  energy  (i.e.,  the  collision  energy)  released 

in  translation.  The  theoretically  calculated  c.m.  angular  distribution  is  very  similar  to  the 

experimentally  derived  angular  distribution  (see  Fig.  17),  but  these  two  distributions  are  not  quite 

in  quantitative  agreement.  The  experimental  result  predicts  slightly  more  forward  scattering  of 
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OC  relative  to  sideways  scattering  than  does  the  theoretical  result.  The  slight  disagreement  in 
the  c.m.  angular  distributions  might  reflect  a  minor  error  in  the  OCO  potential  or,  perhaps,  the 
fact  that  the  calculations  employed  a  classical  rather  than  quantum  scattering  treatment.  The 
good  agreement  between  the  experimental  results  and  the  predictions  of  the  QCT  calculations, 
which  used  only  the  three  low-lying  triplet  potential  energy  surfaces,  suggests  that  intersystem 
crossing  does  not  play  a  role  in  the  dynamics  of  0(3f5)  +  CO  collisions.  The  absence  of 
intersystem  crossing  was  also  a  conclusion  in  the  measurements  of  the  excitation  functions  for 
O (3P)  +  H2  — >  OH  +  H,24  and  O (3P)  +  CH4  — >  OCH3  +  H,9  and  in  the  dynamics  of  the  reaction 

3  8 

o  CP)  +  D2  — >  OD  +D.  The  c.m.  translational  energy  distributions  show  that  reactively 
scattered  OC  tends  to  be  highly  internally  excited.  On  average,  the  reactive  OC  products  are 
formed  with  internal  energies  corresponding  to  roughly  v=5  or  v=6.  But  a  significant  fraction 
(~8  percent)  are  fonned  with  internal  energies  corresponding  to  v=l  1  or  higher.  The  propensity 
of  reactive  OC  products  to  scatter  in  the  forward  direction  might  be  interpreted  as  a 
manifestation  of  a  stripping  reaction,  in  which  the  reaction  occurs  through  high-impact- 
parameter  collisions.  However,  a  typical  characteristic  of  a  stripping  mechanism  is  relatively 
little  energy  transfer  into  product  internal  degrees  of  freedom,  which  is  at  odds  with  the 
observation  of  large  energy  transfers. 

Examination  of  the  calculated  trajectories  provides  insight  into  the  mechanism  that  leads 
to  the  observed  scattering  dynamics.  Several  examples  of  animated  trajectories  may  be  viewed 
online  as  Supporting  Information.  With  the  exception  of  the  large-impact-parameter  collisions 
that  lead  to  inelastic  scattering  with  little  change  in  direction  and  velocity,  the  collisions  that  lead 
to  both  inelastic  and  reactive  products  have  a  striking  similarity.  This  was  shown  earlier  in  the 
extremely  high  degree  of  similarity  of  the  theoretically  calculated  angular  and  translational 
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energy  distributions  for  the  inelastic  and  reactive  channels  without  the  CO(v=0)  contribution. 
The  shallow  wells  of  the  two  lowest  lying  potential  energy  surfaces  facilitate  the  formation  of  a 
transient  OCO  complex,  which  is  formed  over  a  large  range  of  impact  parameters  but  mainly  for 
angles  of  approach  that  are  close  to  the  transition  state  geometry  (O-C-O  angle  of -120°).  The 
complex  can  persist  for  several  vibrational  periods,  but  it  usually  does  not  last  long  enough  for  a 
rotational  period.  Once  the  complex  is  formed,  the  outcome  may  be  an  exchange  of  the  O  atom 
that  is  bound  to  the  C  atom  (reactive  scattering)  or  no  exchange  (inelastic  scattering).  Such  a 
complex  trajectory  that  results  in  inelastic  scattering  may  be  thought  of  as  a  failed  reactive 
collision.  Thus,  for  reactive  scattering  trajectories  and  for  inelastic  scattering  trajectories  that 
transfer  significant  energy  to  the  CO  product,  the  collisions  involve  a  complex  that  tends  to 
channel  CO  products  into  the  forward  direction  but  that  can  also  fling  CO  products  into  the 
sideways  and  backward  directions.  The  details  of  how  the  complex  falls  apart  determine 
whether  or  not  the  outcome  is  reactive  or  inelastic  and  how  the  energy  is  partitioned  into  product 
degrees  of  freedom.  The  theoretical  observation  of  a  transient  OCO  complex  supports  earlier 
suggestions  of  such  a  complex  that  can  be  stabilized  by  collisions,  in  connection  with  isotopic 

25 

fractionation  of  CO2  in  the  atmosphere. 

The  theoretical  calculations  provide  a  detailed  view  of  the  rotational  and  vibrational 
excitation  in  the  CO  products.  Figure  5  shows  that  the  reactive  and  inelastic  rovibrationally 
resolved  cross  sections  extend  to  j — 100.  The  a(0,j)  inelastic  cross  section  peaks  at  much  lower 
rotational  quantum  numbers  and  has  a  much  larger  magnitude,  reflecting  the  different  dynamics 
of  pure  rotational  excitation,  which  occurs  at  large  impact  parameters.  The  maxima  for  the 
inelastic  o(v=l  -3,j)  cross  sections  occur  for  lower  rotational  quantum  numbers  as  the  vibrational 
levels  increase,  reflecting  the  fact  that  less  energy  is  available  for  rotation  as  the  vibrational 
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excitation  increases.  This  trend  is  in  agreement  with  previous  experimental  studies  of  reactive 
triatomic  collisions,-  and  for  experimental  and  theoretical  studies"  ’  on  inelastic  triatomic 
collisions.  In  general,  products  can  be  highly  rotationally  excited,  especially  when  the  product 
vibrational  energy  is  relatively  low. 

As  an  independent  validation  of  the  calculated  rovibrational  distributions,  we  show  in 
Fig.  19  calculated  nascent  emission  spectra  based  on  the  CO(v,j)  populations  for  the  summed 
reactive  and  inelastic  channels  and  measured  nascent  spectra  for  O  +  CO  collisions  near  8  km  s  1 
relative  velocity.  The  calculated  spectra  were  generated  from  the  combined  inelastic  and 
reactive  CO(v,j)  populations  and  tabulated  lifetimes  for  CO(v,j)  states.  Fig.  19a  shows  the 
calculated  results  at  5.0  cm  1  resolution  from  2  to  7  microns,  where  Av=l  (fundamental)  and 
Av=2  (overtone)  bands  are  clearly  visible.  We  note  the  bandhead  structure  in  the  fundamental 
and  overtone  spectra  corresponding  to  R-branch  transitions  from  high  rotational  levels  ( j=50- 
100)  in  excited  vibrational  levels  of  CO.  Within  the  fundamental  band,  each  peak  corresponds  to 
an  R-branch  transition  from  a  particular  vibrational  level,  with  the  shortest  wavelength  peak 
corresponding  to  transitions  from  v=l,  and  so  on.  Fig.  19b  shows  the  same  calculated  nascent 
spectra  degraded  to  0.05  microns  spectral  resolution.  Fig.  19b  also  shows  results  from  two 
measurements,  the  crossed-beams  laboratory  measurements  of  Upschulte  and  Caledonia,11  which 
range  from  4.5  -  5.5  microns,  and  the  space-based  orbital  gas  release  measurements  of  Green  et 
a/.14  (digitized  from  Fig.  12,  17933  scan),  which  range  from  2-7  microns.  The  Green  et  al. 
spectra  capture  the  fundamental  and  overtone  emission.  All  spectra  have  been  normalized  so 
that  their  maxima  are  1.0.  The  spectrum  of  Upschulte  and  Caledonia  was  shifted  by  -0.848 
microns  and  the  spectrum  of  Green  et  al.  was  shifted  by  +0.0772  microns.  There  is  very  good 
agreement  between  all  spectra  for  the  fundamental  band,  except  at  longer  wavelengths  where  the 
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spectrum  of  Green  et  al.  has  apparently  a  low  signal  to  noise  ratio.  The  ratio  of  the  overtone  to 
the  fundamental  peak  heights  for  the  calculated  spectrum  is  about  1 .5  larger  than  the  measured 
spectrum  of  Green  et  al.  This  difference  may  be  the  result  of  uncertainties  in  combining  spectra 
from  two  different  circular  variable  filters  (CVFs)  covering  different  wavelength  regions  (0.7-5 
microns  and  4.7-22  microns),  other  uncertainties  involving  the  absolute  determination  of  relative 
velocities  of  O  and  CO,  the  sensitivity  of  the  measured  spectrum  to  the  direction  of  the  CVF  field 
of  view,  other  factors  in  these  difficult  space  measurements,  or  in  deficiencies  in  the  present 
calculations.  More  highly  resolved  spectra  are  required  to  detennine  the  origins  of  these 
differences. 

The  current  experimental  and  theoretical  results  reveal  similar  limits  in  the  amount  of 
internal  excitation  in  scattered  CO  to  what  was  reported  in  the  IR  emission  experiment  conducted 
by  Upschulte  and  Caledonia.  They  concluded  that  collisions  of  0(  P)  with  CO  at  Eco\\  =  77.9 
kcal  mol  1  can  produce  CO  in  vibrational  levels  up  to  v=ll,  while  our  crossed-beams  study  of 
Of3/3)  with  ClsO  at  (Ccoii)  =  83  kcal  mol'1  verified  that  both  the  inelastic  product,  C180,  and  the 
reactive  product,  16OC,  may  be  in  vibrational  levels  of  v=ll  or  greater.  A  high  degree  of 
consistency  is  emerging  between  the  present  measurements,  the  spectral  measurements  of 
Caledonia  and  Upschulte,  the  spectral  measurements  of  Green  et  al.,  and  the  theoretical 
calculations,  especially  with  regard  to  the  high  degree  of  internal  excitation  of  the  product  CO. 

V.  CONCLUDING  REMARKS 

Crossed-molecular  beams  methods  and  theoretical  calculations  have  been  used  to 
investigate  the  hyperthermal  interactions  of  Of  P)  with  CO  at  c.m.  collision  energies  near  83 
kcal  mol  1  (corresponding  to  a  relative  velocity  near  8  km  s_1).  Inelastic  scattering  experiments 
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on  the  process,  160(3P)  +  C180(v,j)  — >  C180(v',j')  +  160,  indicate  that  oxygen  atoms  mainly 
scatter  from  C  O  with  relatively  large  impact  parameters  and  lead  to  forward  scattering  with 
little  change  in  direction  and  translational  energy.  Averaged  over  all  inelastic  collisions,  ~15 
percent  of  the  collision  energy  is  transferred  to  internal  energy  in  CO.  Although  this  fraction  is 
relatively  small,  the  collision  energy  is  high,  so  internal  excitation  of  CO  above  an  energy 
corresponding  to  v=l  still  occurs  in  approximately  84  percent  of  all  inelastic  collisions.  The 
reactive  channel,  160 (3P)  +  C180— >■  16OC  +  180,  created  16OC  products  that  are  frequently 
scattered  in  the  forward  direction,  although  the  fraction  of  collisions  that  lead  to  16OC  scattered 
in  the  forward  direction  is  comparable  to  the  fraction  that  lead  to  sideways  and  backward 
scattering.  The  16OC  product  formed  in  the  reactive  channel  had  a  very  broad  translational 
energy  distribution,  with  ~40  percent,  on  average,  of  the  available  (collision)  energy  transferred 
into  internal  excitation.  Approximately  97  percent  of  reactive  collisions  produce  16CO  with 
internal  energies  greater  than  an  energy  corresponding  to  v=l,  and  about  8  percent  of  collisions 
lead  to  internal  energies  greater  than  v=ll.  Quasiclassical  trajectory  calculations  on  computed 
triplet  potential  energy  surfaces  are  in  excellent  agreement  with  the  experimental  results.  The 
calculations  reveal  a  mechanism  for  large  energy  transfers,  into  both  rotation  and  vibration, 
where  the  reagent  O  atom  encounters  the  carbon  end  of  the  CO  molecule  a  range  of  impact 
parameters  in  a  bent  O-C-O  geometry  (~80°-140°).  The  three  atoms  linger  together  for  less  than 
a  rotational  period  but  usually  for  more  than  a  CO  vibrational  period,  and  the  transient  complex 
breaks  up  into  either  inelastic  or  reactive  products.  The  inelastic  trajectories  that  lead  to 
significant  energy  transfer  to  the  CO  may  thus  be  considered  failed  reactive  collisions.  The 
similarity  in  inelastic  and  reactive  trajectories  that  lead  to  high  vibrational  excitation  in  the  CO 
product  is  supported  by  calculated  translational  energy  and  angular  distributions  for  reactive  and 
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inelastic  collisions  without  the  CO(v=0)  product.  These  distributions,  which  remove 
contributions  from  purely  inelastic  rotational  excitation,  are  nearly  identical.  The  theoretical 
calculations  have  been  used  to  predict  inelastic  and  reactive  cross  sections  at  a  variety  of 
collision  energies  and  to  predict  infrared  emission  spectra  that  compares  favorably  with  earlier 
measurements.  Given  the  accuracy  with  which  the  calculations  predict  the  collision  dynamics 
observed  in  our  high-fidelity  crossed-beams  experiments,  we  expect  the  mechanistic  insights  and 
additional  predictions  offered  by  the  theory  to  be  accurate,  too. 
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SUPPORTING  INFORMATION 

Representative  animated  trajectories  of  O  +  CO  collisions  are  available  as  Supporting 
Information.  These  trajectories  represent  exchange  of  the  oxygen  atom  that  is  bound  to  the 
carbon  atom  (reactive  scattering)  and  collisional  excitation  of  CO  without  exchange  of  the 
oxygen  atom  (inelastic  scattering).  Some  trajectories  show  the  formation  of  a  transient  complex 
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in  which  the  three  atoms  linger  together  for  a  short  time,  usually  less  than  a  rotational  period, 
before  separating  into  products.  This  material  is  available  free  of  charge  via  the  Internet  as  a 
Link  associated  with  the  manuscript. 
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FIGURE  CAPTIONS 


Figure  1.  Calculated  potential  energy  curves  of  the  ground  singlet  state  of  CO2,  ( lA'),  and  the 
three  lowest  triplet  states, 3 A',  1  3 A",  and  2  A",  at  an  O-C-O  angle  of  120°  and  where  one  of  the 
C-0  bond  distances  is  fixed  at  1 .2  A  and  the  other  is  varied.  The  curves  show  these  states,  at  Cs 
geometries  to  their  O  +  CO  dissociation  limits,  as  a  function  of  the  variable  C-0  distance,  Rco- 

Figure  2,  Calculated  differential  energy  cross  sections,  do/dEtmns,  for  O  +  CO  collisions  at  8  km 
s  '  collision  velocity  (Ec 0n  =  77.9  kcal  mol-1).  —  (black  line)  total  reactive  channel,  —  (red  line) 
total  inelastic  channel,  —  (black  dotted  line)  total  reactive  minus  contribution  from  reactive 
CO(v=0)  product,  —  (red  dotted  line)  total  reactive  minus  contribution  from  reactive  CO(v=0) 
product. 


Figure  3.  Calculated  differential  angular  cross  sections,  do/dtl,  for  O  +  CO  collisions  in  the 
c.m.  frame  at  8  km  s_1  collision  velocity  (Eco\\  =  77.9  kcal  mol  ’).  —  (black  line)  total  reactive 
channel;  —  (red  line)  total  inelastic  channel,  —  (black  dotted  line)  total  reactive  minus 
contribution  from  reactive  CO(v=0)  product,  —  (red  dotted  line),  total  reactive  minus 
contribution  from  reactive  CO(v=0)  product.  For  the  reactive  channel,  the  c.m.  scattering  angle, 
©cm,  is  defined  as  the  inverse  cosine  of  the  dot  product  of  the  unit  velocity  vectors  in  the  c.m. 
frame  of  the  incoming  (reagent)  O-atom  and  outgoing  (product)  CO  molecule.  For  the  inelastic 
channel,  the  c.m.  scattering  angle,  0cm,  is  defined  as  the  inverse  cosine  of  the  dot  product  of  the 
unit  velocity  vectors  in  the  c.m.  frame  of  the  incoming  (reagent)  O  atom  and  outgoing  (product) 
O  atom. 


Figure  4.  Calculated  cross  sections  for  O  +  CO  collisions  at  8  km  s^1  collision  velocity  (Eco\\  = 
77.9  kcal  mol-1),  showing  the  contribution  from  each  final  CO(v,j)  state  as  a  function  of  the  final 
translational  energy  of  collision  products,  (a)  reactive  channel,  (b)  inelastic  channel,  (c)  inelastic 
channel  minus  the  CO(v=0,j)  product  contributions.  The  numbers  in  the  figures  indicate  the 
thresholds  for  the  CO(v)  product  vibrational  levels.  In  (b),  the  contribution  from  CO(v=0)  below 
the  CO(v=l)  threshold  extends  above  the  axis  limit  and  so  is  not  visible  in  the  figure. 


Figure  5.  Calculated  rovibrationally  resolved  cross  sections  for  O  +  CO  collisions  at  8  km  s_1 
collision  velocity  (Ccoii  =  77.9  kcal  moF1).  (a)  reactive  channel,  (b)  inelastic  channel.  —  (black 
line)  CO(v=0),  —  (red  line)  CO(v=l),  —  (blue  line)  CO(v=2),  —  (green  line)  CO(v=3).  The 

_ o 

CO(v=0)  inelastic  cross  section  has  been  multiplied  by  1  x  1 0  in  order  to  fit  on  the  scale. 
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Figure  6.  Calculated  vibrationally  resolved  cross  sections  for  O  +  CO  collisions  at  3,  4,  5,  6,  7, 
8,  9,  10  km  s-1  collision  velocities,  (a)  reactive  collisions,  (b)  inelastic  collisions.  The  numbers 
in  the  figure  refer  to  the  collision  velocity. 


Figure  7.  Calculated  CO  product  energy  fractions  for  O  +  CO  collisions  as  a  function  of 
collision  energy,  o  reactive  channel,  □  inelastic  channel,  —  vibrational  energy  fraction  (f v), 
—  rotational  energy  fraction  (fJ). 


Figure  8.  Calculated  vibrational  excitation  and  photon  cross  sections  for  O  +  CO  collisions  as  a 
function  of  collision  energy,  o  reactive  channel,  □  inelastic  channel,  —  results  for  the  photon 

cross  section  ap  =  a(v=l)  +  2a(v=2)  +  3a(v=3)+..., - results  for  the  vibrational  excitation 

function  ae  =  a(v=l)  +  a(v=2)  +  a(v=3)  +  . . .  . 

Figure  9.  Pictorial  diagram  of  the  crossed  molecular  beams  apparatus. 

Figure  10.  Translational  energy  distributions  of  atomic  and  molecular  oxygen  in  the 
hyperthennal  beam.  The  mole  fraction  of  O  atoms  is  98  percent,  and  the  mole  fraction  of  O2 
molecules  is  2  percent. 

Figure  11.  Newton  diagram  for  collisions  of  0(  P)  with  C  O  at  £con  =  83  kcal  mol  ,  where  O 
has  a  velocity  of  8100  m  s_1  and  the  C180  beam  velocity  is  800  m  s_1.  The  shaded  region 
indicates  the  range  of  laboratory  angles  examined.  The  radii  of  the  blue  circles  represent  the 
recoil  velocities  of  elastically  scattered  C180  (solid)  and  160  (dashed).  The  radii  of  the  red 
circles  represent  the  maximum  recoil  velocities  of  reactively  scattered  16OC  (solid)  and  180 
(dashed). 

Figure  12.  Time-of-flight  distributions  of  inelastically  scattered  O  and  C  O  following 
collisions  with  (Ec 0n)  =  83  kcal  mol  '.  The  circles  are  the  experimental  data.  The  red  (O)  and 
blue  (C  O)  solid  curves  are  the  forward-convolution  simulations  of  the  data,  derived  from  the 
c.m.  angular  and  translational  energy  distributions  of  corresponding  colors  in  Fig.  14. 


1 8 

Figure  13.  Laboratory  angular  distributions  of  O  atoms  and  C  O  molecules  that  scattered 
inelastically  with  (Ec 0n)  =  83  kcal  moT1.  The  circles  with  error  bars  are  the  experimental  data, 
and  the  colored  lines  are  the  forward-convolution  fits  to  the  data,  derived  from  the  c.m.  angular 
and  translational  energy  distributions  of  corresponding  colors  in  Fig.  14.  The  error  bars  are 
estimated  from  fitting  the  experimental  TOF  distributions  with  a  modified  Gaussian  function  and 
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finding  areas  of  the  maximum  and  minimum  acceptable  fits  by  adjusting  the  Gaussian 
parameters.  The  error  bars  thus  represent  the  maximum  and  minimum  integrals  of  the  TOF 
distributions  based  on  our  best  judgment  and  are  expected  to  be  greater  than  ±2c. 

Figure  14.  Center-of-mass  angular  and  translational  energy  distributions  for  the  inelastic 
scattering  of  O  with  C  O  at  (Eco\\)  =  83  kcal  mol  ,  used  in  the  forward-convolution  fit  of  the 
laboratory  TOF  and  angular  distributions  for  m/z  =16  (0+)  and  30  (C180+).  These  distributions 
were  obtained  from  QCT  calculations  for  O  (red)  and  C  O  (blue). 

Figure  15.  Time-of-flight  distributions  of  reactively  scattered  OC  following  reaction  of  Of  P) 
with  C  O  at  (ficoii)  =  83  kcal  mol  .  The  circles  are  the  experimental  data.  The  dashed  blue  and 
solid  red  curves  are  the  forward-convolution  simulations  of  the  data,  derived  from  the  c.m. 
angular  and  translational  energy  distributions  of  corresponding  colors  in  Fig.  17.  The  dashed 
blue  curves  are  predicted  by  theory,  and  the  solid  red  curves  are  the  results  from  optimizing  the 
c.m.  distributions  (shown  in  Fig.  17)  to  fit  the  data.  The  experimental  (solid  red)  and  the 
theoretical  (dashed  blue)  curves  are  nearly  identical  for  each  TOF  distribution. 

Figure  16.  Laboratory  angular  distribution  (circles  with  error  bars)  of  the  reactively  scattered 
16OC  product  following  reaction  of  O (3P)  with  C180  at  (Ecou)  =  83  kcal  mol-1.  The  dashed  blue 
and  solid  red  curves  are  the  forward-convolution  simulations  of  the  data,  derived  from  the  c.m. 
angular  and  translational  energy  distributions  of  corresponding  colors  in  Fig.  17.  The  dashed 
blue  curve  is  predicted  by  theory,  and  the  solid  red  curve  is  the  result  from  optimizing  the  c.m. 
distributions  (shown  in  Fig.  17)  to  fit  the  data.  The  error  bars  are  estimated  from  fitting  the 
experimental  TOF  distributions  with  a  modified  Gaussian  function  and  finding  areas  of  the 
maximum  and  minimum  acceptable  fits  by  adjusting  the  Gaussian  parameters.  The  error  bars 
thus  represent  the  maximum  and  minimum  integrals  of  the  TOF  distributions  based  on  our  best 
judgment  and  are  expected  to  be  greater  than  ±2c. 

Figure  17.  Center-of-mass  angular  and  translational  energy  distributions  for  reactive  scattering 
of  Of’/3)  with  C180  at  (Eco\\)  =  83  kcal  mol-1,  used  in  the  forward-convolution  fit  of  the 
laboratory  TOF  and  angular  distributions  for  m/z  =  28  (16OC+).  The  angular  distribution  (left) 
pertains  to  the  reactively  scattered  16CO.  The  distributions  shown  in  dashed  blue  were  obtained 
from  QCT  calculations,  and  those  shown  in  red  were  determined  by  optimizing  the  c.m. 
distributions  to  fit  the  data. 

Figure  18.  (a)  Experimental  and  (b)  theoretical  c.m.  velocity-flux  maps  for  the  reactive  16OC 
product,  derived  from  the  angular  and  translational  energy  distributions  shown  in  Fig.  17. 

Figure  19.  Infrared  emission  spectra  from  nascent  CO  after  hyperthermal  O  +  CO  collisions, 
(a)  Calculated  spectrum  at  5  cm  1  spectral  resolution  for  collisions  at  8  km  s^1  collision  velocity 
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(£coii  =  77.9  kcal  mol  ').  (b)  —  calculated  spectrum  as  in  (a)  degraded  to  0.05  pm  spectral 

resolution,  measured  spectrum  from  Green  et  al6, - measured  spectrum  of  Upschulte  and 

Caledonia.5  The  spectrum  of  Green  et  al6  has  been  shifted  by  +0.0772  pm,  and  the  spectrum  of 
Upschulte  and  Caledonia5  has  been  shifted  by  -0.0848  pm.  All  spectra  are  normalized  so  that 
the  maximum  intensity  value  is  1.0. 
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SUPPORTING  INFORMATION 


Representative  animated  trajectories  of  O  +  CO  collisions  are  available  as  Supporting  Information. 
These  trajectories  represent  exchange  of  the  oxygen  atom  that  is  bound  to  the  carbon  atom  (reactive 
scattering)  and  collisional  excitation  of  CO  without  exchange  of  the  oxygen  atom  (inelastic  scattering  or 
“non-reactive”  scattering).  Some  trajectories  show  the  formation  of  a  transient  complex  in  which  the 
three  atoms  linger  together  for  a  short  time,  usually  less  than  a  rotational  period,  before  separating  into 
products.  This  material  will  be  available  free  of  charge  via  the  Internet  as  a  Link  associated  with  the 
manuscript. 


General  description: 

0(3P)  +  CO  collision  trajectory  movies  at  8  km  s'1: 

The  molecular  potential  and  methods  are  the  same  as  M.  Braunstein  and  J.  Duff, 

"Electronic  Structure  and  Dynamics  of  O  +  CO  Collisions",  J.  Chem.  Phys.,  1 12,  2736-2745  (2000). 

The  format  is  .avi.  Each  file  contains  ~10  trajectories  and  takes  about  two  minutes  of  wall  clock  minutes 
to  view.  There  are  30  files  total  in  two  folders  for  a  total  of  ~560  Mbytes  (uncompressed). 


Detailed  description  and  notation  of  movie  plots: 

O  +  CO  Collision  velocity  8  km  s'1 

reactive  all  folder  ->  all  reactive  collisions  regardless  of  final  CO(vibrational  level 
nonreactivevgtO  folder  ->  non  reactive  (inelastic)  collisions  for  CO(v>0)  products 

red  sphere:  incoming  O-atom 

black  sphere:  carbon  atom 

blue  sphere:  oxygen  atom  in  target  CO 

0*(incoming)  ->  C— O(target) 

red  line  0*C  distance 
green  line  C— O  distance 
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blue  line  O* — O  distance 


bl->  3 A'  state  (lowest  OCO  triplet  state) 

al  1 3.A  state  (second  lowest  OCO  triplet  state) 

Each  movie  frame  (sphere  symbol)  takes  4.8e-15  seconds  (4.8  fs  =  200  au  time). 
Each  trajectory  is  followed  for  approximately  1.25e-13  seconds  (125  fs) 

The  vibrational  period  of  CO  ranges  from  ~10  fs  (CO(v=l))  to  ~1  fs  (CO(v=10)) 
The  rotational  period  of  CO  ranges  from  ~l.e4  fs  (CO(j=l))  to  —100  fs  (C0(j=100)) 

All  collisions  occur  at  8  km  s-1  center  of  mass  relative  velocity  (77.9  kcal  mol-1) 
0(3P)  +  C0(v=0,j=300  K). 

complex  time  =>  when  rl  +  r2  +  r3  (red  line,  green  line,  blue  line)  <12  a.u. 

fs  =  femtoseconds  (l.e-15  seconds) 

angle  (center  of  mass  scattering  angle) 

v#  =  final  vibrational  state 

j#  =  final  rotational  state 

etrans  =  final  translational  energy  (kcal  mol-1) 

Representative  movie  frames  are  shown  below. 
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Figure  1.  A1  Reactive  trajectory 
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Trajectory  #2 
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Figure  2.  A1  Nonreactive  (inelastic)  trajectory 
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Figure  3.  B1  Reactive  trajectory 
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Trajectory  #2 
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Figure  4.  B1  Nonreactive  (inelastic)  trajectory 
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